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ABSTRACT 

Aims. The condensation of diffuse gas into molecular clouds and dense cores occurs at a rate driven largely by turbulent 
dissipation. This process still has to be caught in action and characterized. 

Methods. We observed a mosaic of 13 fields with the IRAM-PdB interferometer (PdBI) to search for small-scale structure 
m the ^2 CO (1-0) line emission of the turbulent and translucent environment of a low-mass dense core in the Polaris 
Flare. The large size of the mosaic (l'x2') compared to the resolution (4") is unprecedented in the study of the small- 
scale structure of diffuse molecular gas. 

Results. The interferometer data uncover eight weak and elongated structures with thicknesses as small as ~ 3 mpc 
(600 AU) and lengths up to 70 mpc, close to the size of the mosaic. These are not filaments because once merged 
with short-spacings data, the PdBl-structures appear to be the sharp edges, in space and velocity-space, of larger-scale 
structures. Six out of eight form quasi-parallel pairs at different velocities and different position angles. This cannot 
be the result of chance alignment. The velocity-shears estimated for the three pairs include the highest values ever 
measured in regions that do not form stars (up to 780 km s~^ pc~^). The CO column density of the PdBTstructures 
is in the range A'^(CO) = 10^'* to 10^^ cm~^ and their H2 density, estimated in several ways, does not exceed a few 10"^ 
cm~^. Because the larger scale structures have sharp edges (with little or no overlap for those that are pairs), they have 
to be thin layers of CO emission. We call them SEE(D)S for Sharp-Edged Extended (Double) Structures. These edges 
mark a transition, on the milliparsec scale, between a CO-rich component and a gas undetected in the ^^CO(l-O) line 
because of its low CO abundance, presumably the cold neutral medium. 

Conclusions. We propose that these SEE(D)S are the first directly-detected manifestations of the intermittency of 
interstellar turbulence. The large velocity-shears reveal an intense straining field, responsible for a local dissipation rate 
several orders of magnitude above average, possibly at the origin of the thin CO layers. 

Key words. ISM: evolution - ISM: kinematics and dynamics - ISM: molecules - ISM: structure - ISM: general - Turbulence 
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1. Introduction 

Turbulence in the interstellar medium (ISM) remains a puz- 
zle in spite of dedicated efforts on observational and nu- 
merical grounds. This is because it is compressible, magne- 
tized, and multi-phase, but also because of the huge range 
of scales separating those of injection and dissipation of 
energy. Moreover, because turbulence and magnetic fields 
are the main support of molecular clouds against their self- 
gravity, turbulent dissipation is a key process among all 
th ose eventually leading to star formation (see the re views 
of lElmegreen fc Scald[200l iScalo fc Elmegreenl[200l . 

In molecular clouds, turbulence is observed to be highly 
supersonic with respect to the cold gas. It is thus antici- 
pated to dissipate in shocks in a cloud-crossing time (i.e. ~ 
a few 10 Myr for giant molecular clouds of 100 pc with inter- 
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Fig. 1. The location of the 13-field mosaic observed 
at the Plateau de Bure interferometer (centered at 
RA=01:55:12.26 and Dec=87:41:56.30) is shown as the box 
on top of the integrated emission of the ^^CO and "'^'^CO 
(J=l-0) maps obtained at the lRAM-30m. This is a place 
of low, almost featureless, CO line brightness. The arc-like 
structure visible in ^^CO traces the outer layers of the low- 
mass dense core. Contour levels are shown in the wedges. 
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nal velocity dispersion of a few km s~^). M agnetic fields do 
not s ignificantly slow the dissipation down ()Mac Low et al.l 
119981) . Actually this is the basis of the tur bulent models of 
star formation (jMac Low fc Klessenll200j ) - one of the two 
current scenarii of low-mass star formation - in which self- 
gravitating entities form in the shock-compressed layers of 
supersonic turbulence. 

However, while it is unquestionable that the ISM 
is regularly swept by large-scale shock-waves triggered 
by supernovae explosions th at partly feed t h e in- 
terstellar turbulent cascade (iJoung fc Mac Low! 120061 
Ide Avillez fc Breitschwerdtl l2007| ). the smallest scales, 
barely subparsec in these simulations, are still orders of 
magnitude above the smallest observed structures and are 
unlikely to provide a proper description of the actual 
dissipation processes. Whether turbulent dissipation oc- 
curs primarily in compressive (curl-free) or in solenoidal 
(divergence-free) modes in the interstellar medium has 
therefore to be considered as an open issue. 

An ideal target to study turbulent dissipation is the 
diffuse molecular gas because it is the component in which 
dense cores form, with less turbulent energy density than 
their environment. The word "diffuse" here comprises all 
material in the neutral ISM at large that is not in dense 
cores i.e. whose total hydrogen column density is less than a 
few 10^^ cm^^. This includes the mixture of cold and warm 
neutral medium (CNM and WNM), the edges of molecular 
cloud complexes (also called translucent gas), and the high 
latitude clouds. Diffuse gas builds up a ma jor mass fraction 
of th e ISM. Actually, on the 30 pc scale, I Goldsmith et al.l 
(|2008i)find that half the mass of the Taurus- Auriga-Perseus 
complex lies in regions having H2 column density below 
2.1 X 10^1 cm-2. 

Turbulent dissipation may also provide clues to the 
"outstanding mysteries" raised by observations of d iffuse 
molecular gas (see the review o f[S now fc McCal]|[200l: the 
ubiqu itous small scale structure, down to AU-scales (|Heilesl 
|2007() . the remarkable molecular richness found in this 
hostile medium, weakly shielded from UV radiation (e.g. 
iLiszt fc LucaslllQQSHGredel et al.ll2002[ ). the bright emission 
in the H2 pure rotational lines exceeding t he predictions of 
photon-dominated region (PDR) models ([Falgarone et al.l 
l2005t iLacour et a"L| [2005). the ^^CO small-scale structures 
with a broad range of temperatures, H2 densities and 
linewidths that precl ude a single interpretation in terms 
of co l d deiise clumps flngalls et al."200oi 120071 : iHeithausenI 
12004 1200d fSakamoto fc Smiada 2003). 

The present paper extends the investigation of of tur- 
bulence down to the mpc-scale in the translucent environ- 
ment of a low-mass dense core of the Polaris Flare. Over the 
years, this investigation has progressed along three comple- 
mentary directions: 

(i) A two-point statistical analysis of the velocity 
field traced by the ^^CO line emission, and conducted 
on maps of increasing size. Using n umerical s imula - 
tion s of mild ly compressible turbulence, iLis et all ()1996[ ) 
and iPetv fc Falgarone ( 2003) first proposed that the non- 
Gaussian probability distribution functions (pdjs) of line 
centroid velocity increments (CVI) be the signatures of the 
space-time intermittency of turbulence Q because the ex- 



trema of CVI (E-CVI) trace extrema of the line-of-sight 
average of the modulus of the pos vorticity. Statistical 
analysis conducted on parsec-scale maps in two nearby 
molecular clouds have revealed t hat these extrema form 
parsec-scale coherent struct ures (|Hilv-Blant et al.l 120081 : 
iHilv-Blant fc Falgaronel[2009l resp. Paper III, HF09). 
(ii) A detailed analysis (density, temperature, molecular 
abundances) of these coherent structures, based on their 
molecular line emission. The gas there is more optically thin 
in the ^^C O lines, warmer and more dilu te than the bulk 
of the gas (jHilv-Blant fc Falgarone! [20071 hereafter Paper 
II), and large HCO"*" abundances, unexpected in an envi- 
ronment wea kly shielded from UV radiation, have been de- 
tected there ([Falgarone et al.[[2006l Paper I). 
(Hi) Chemical models of non-equilibrium w arm chemistry 
trigg ered by bursts of turbulent dissipation ([Joulain et al.l 
Il998[ ). The most recent progresses along those lines in- 
clude the c hemical model s of tu rbulent dissipation regions 
(TDRs) bv iGodard et al.l ([2009( ) and their successful com- 
parison to several data sets, among which new submillime- 
ter detections of ^^CII+(l-0) (Falgarone et al. in prepara- 
tion). 

The ^^C0(J=2-1) observations of the Polaris Flare with 
unprecedented angular resolution and dynamic range are 
the first to evidence the association between extrema of 
CVI and observed velocity-shearj[3 (HF09). No shock sig- 
nature (density and/or temperature enhancement, SiO de- 
tection) has been found in the coherent structure of E-CVI 
identified in the Polaris Flare (Hily-Blant and Falgarone, in 
preparation). All the above suggest (but docs not prove yet) 
that the coherent structures carrying the statistical prop- 
erties of intermittency are regions of intense velocity-shears 
where dissipation of turbulence is concentrated. 

The CO (1-0) observations reported in this paper have 
been performed in a field located on one branch of the 
Polaris Flare E-CVI structure, in the translucent and fea- 
tureless environment of a dense core (Fig.[T]). The outline of 
the paper is the following: the observations and data reduc- 
tion are described in Section 2. The observational results 
are given in Section 3. The characterization of the emit- 
ting gas is made in Section 4 and we discuss, in Section 5, 
the possible origin and nature of the CO structures that we 
have discovered. Section 6 puts our results in the broad per- 
spective provided by other data sets and Section 7 compares 
them to chemical model predictions and numerical simula- 
tions of turbulence. The conclusions are given in Section 



Intermittency here refers to the empirical property of 
high Reynolds number turbulence to present an excess of rare 
events compared to Gaussian statistics, this excess being in- 



creasingly large as velocity fluctuations at smaller and smalle r 
scales are considered (see the review of lAnselmet et al.|[200I^ . 
Alt hough the origin of intermittency is st i ll an open issue (but 
see iMordant et al'l 120021 : iChevillard et all 120051 : lArneodo etlll 
l2ooa r it is quantitatively characterized by the anomalous scal- 
ing of the high-order structure functions of the velocity and 
the shape of non-Gaussian pdjs of quantities involving velocity 
derivatives (e.s. iFrisch 1995). 

^ We use velocity-shear rather than velocity-gradient because 
the observations provide cross-derivatives of the velocity field, 
i.e. the displacement measured in the plane-of-the-sky (pos) is 
perpendicular to the line-of-sight velocity 
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2. IRAM Plateau de Bure Interferometer 
observations 

We used the IRAM Plateau de Bure Interferometer (PdBI) 
to image, at high angular resolution and in the ^^CO (J=l- 
0) line, a region of ^ 1' x 2' in the tra nslucent environ- 
ment of a dense core in the Polaris Flare (|Heithausenlll999t 
iHeithausen et all l2002t ). The location of the target field 
is shown in Fig. [T] as a rectangle on larger scale, single- 
dish maps of i ntegrated ^^CO(J=1- 0) and 13CO(J=1-0) 
emission from ()Falgarone et all 119981 hereafter F98). The 
average column density in this region (~ 10^^ cm~^) is 
about 100 times smaller than in the central parts of the 
dense core (~ 10^"^ cm~^), 3 arcmin westwards. The aver- 
age integrated ^'^CO intensity over the mosaic area is weak 
W(^3co) = 2 Kkms"^ 

2.1. Observations 

The observations dedicated to this project were carried 
out in 1998 and 1999 with the IRAM interferometer at 
Plateau de Bure in the C and D configurations (base- 
line lengths from 24m to 161m). One correlator band of 
10 MHz was centered on the ^^CO (J=l-0) frequency to 
cover a ~ 23kms~^ bandwidth with a channel spacing of 
39 kHz, i.e. O.lkms"^. Four additional correlator bands 
of 160 MHz were used to measure the 2.6 mm continuum 
over the 500 MHz instantaneous IF-bandwidth then avail- 
able. 

We observed a 13-field mosaic centered on q;2ooo = 
0l''55"12.26^ (52000 = 87°41'56.30". The field po^sitions 
followed a compact hexagonal pattern to ensure Nyquist 
sampling in all directions and an almost uniform noise over 
a large fraction of the mosaic area (see Fig. lA.ll of Appendix 
A). The imaged field-of-view is about a rectangle of dimen- 
sion of 85" X 130" oriented at a position-angle of 15° (be- 
cause the (RA,Dec) PdBI field was selected in maps made 
in (l,b) coordinates). 

Polaris being a circumpolar source, this project was a 
good time-filler. It was thus observed at 22 different occa- 
sions, giving a total of about 180 hours of telescope time 
with most often 3 or 4 antennas and rarely 5 antennas. 
Taking into account the time for calibration and data filter- 
ing this translates into on-source integration time of useful 
data of 65.2 hours for a full 5-antenna array. The typical 
2.6 mm resolution of these data is 4.3". The data used to 
produce the missing short-spacings are those of the IRAM 
key-program, fully described in F98 (see also Table 1). 

2.2. Data reduction 

The data processing w as done with the GILDAS0 soft- 
ware suite (|Petvll2005l ). Standard calibration methods im- 
plemented in the GILDAS/CLIC program were apphed us- 
ing close bright quasars as calibrators. The calibrated uv 
tables were processed through an Banning filter which spec- 
trally smoothed the data (to increase the intensity signal- 
to-noise ratio) while keeping the same velocity /frequency 
channel spacing. 

All other processing took place into the 
GILDAS/MAPPING software. Following iGueth et all 

^ See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS softwares. 



((l99l . the single-dish map from the IRAM-30m key pro- 
gram were used to create the short-spacing visibilities not 
sampled at the Plateau de Bure. These were then merged 
with the interferometric observations. Two different sets of 
uv tables (i.e. with and without short-spacings) were then 
imaged using the same method. Each mosaic field was 
imaged and a dirty mosaic was built combining those fields 
in th e following o ptimal way in terms of signal-to-noise 
ratio (|Guethll200ll) 

J{a, 5) = ^2 ^»(«' ^) I 2^, ^2 • 

In this equation, J(a, (5) is the brightness distribution in the 
dirty mosaic image, Bi are the response functions of the i 
primary antenna beams, Fi are the brightness distributions 
of the individual dirty maps and Oi are the corresponding 
noise values. As may be seen in this expression, the dirty 
intensity distribution is corrected for primary beam attenu- 
ation, which makes the noise level spatially heterogeneous. 
In particular, noise strongly increases near the edges of the 
field of view. To limit this effect, both the primary beams 
used in the above formula and the resulting dirty mosaics 
are truncated. The standard level of truncation is set at 
20% of the maximum in MAPPING. 

Deconvolution methods were different for both data sets 
(i.e. with and without short-spacings). The dirty image of 
the PdBI-only data was deconvolved using the standard 
Clark CLEAN algorithm. One spatial support per channel 
map was defined by selecting positive regions on the first 
clean image which was obtained without any constraint. 
This geometrical constraint was then used in a second de- 
convolution. While it can bias the result, this two-step pro- 
cess is needed when deconvolving interferometric observa- 
tions of extended sources without short-spacings. Indeed, 
the lack of short-spacings implies (among other things) a 
zero valued integral of the dirty beam and dirty image, 
which in turn perturbs the CLEAN convergence when the 
source is extended because the algorithm searches as much 
positive as negative CLEAN components. The only way 
around is to guide the deconvolution by the definition of a 
support where the signal is detected. On the other hand, 
the deconvolution of the combined short-spacings and in- 
terferometric uv visibilities can process blindly (i.e. without 
the possible bias of defining a support where to search for 
CLEAN components). This is what has been done and the 
good correlation of the structures seen in the deconvolved 
images of the data with and without short-spacings (see 
Fig. 4 and 5) gives us confidence in our deconvolution of 
the PdBI-only data. 

The two resulting data cubes (with and without short- 
spacings) were then scaled from Jy/beam to T„ib temper- 
ature scale using the synthesized beam size (see Table [1]). 
Final noise rms measured at the centered of the mosaic is 
about 0.23 K in both data cubes. 



3. Observational results 

3.1. PdBI structures: sharp edges of extended structures 

At the adopted cloud distance of d — 150 pc, 1 " corre- 
sponds to 0.75 mpc or 150 AU, so that the spatial resolution 
of the PdBI data is 3.2 mpc or 660 AU. 
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Fig. 3. From top to bottom, maps of the PdBI, PdBI+SOm and 30m of ^^CO(l-O) emission integrated over the same 
velocity shces of 0.3 km centered as indicated. 
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Table 1. Observation parameters. The projection center of all the data displayed in this paper is: 02000 = 0l''55™12.26^, 
^2000 = 87°41'56.30". 



Molecule 


Transition 


Frequency 


Instrument 


Config. 


Beam 


PA 


Vel. Resol. 


Int. Time 


Noise 






GHz 






arcsec 





km s"'^ 


hours 


K 




(J=l-0) 


115.271195 


PdBI 


C&D 


4.4 X 4.2 


80 


0.1 


65.2/180" 


0.23" 






115.271195 


30m 




21.3 





0.1 


-/- 


0.40 


13CO 


(J=l-0) 


110.201354 


30m 




22.3 





0.1 


-/- 


0.19 



" Two values are given for the integration time: the 5 antennae array equivalent on-source time and the telescope time. 
' The noise value quoted here is the noise at the mosaic phase center 



Table 2. Spatial and kinematic characteristics of the ^^CO PdBI-only structures. 



Structure 


km s^^ 


km s~^ 


Aui/2 

km s^""^ 


K 


W{CO) 
Kkms-i 


dl/2 

arcsec 


«± 

mpc 


PA 




max b 

cm 


1 


-5.7 


-5.2 


0.1 


0.6 


0.06 


4 


3.0 


109 


1000 


2 


-5.6 


-5.4 


0.2 


1.8 


0.36 


10 


7.5 


173 


2400 


3 


-5.2 


-4.8 


0.2 


2.4 


0.48 


9 


6.8 


62 


3200 


4 


-4.3 


-4.1 


0.1 


1.2 


0.12 


8 


6.0 


59 


1000 


5 


-3.4 


-2.6 


0.4 


4 


1.6 


12 


9.0 


91 


8900 


6 


-3.4 


-2.6 


0.25 


1.2 


0.3 


10 


7.5 


161 


2000 


7 


-3.2 


-3.0 


0.15 


1.2 


0.18 


15 


11.3 


173 


800 


8 


-1.7 


-1.3 


0.15 


1.2 


0.18 


9 


6.8 


59 


1200 



° projected thickness of the filamentary structures deconvolved from beam size 

^ upper limit because computed as riHa = N{]i2)/li_ instead of using with Ai'(H2) derived from W(CO) (see text) 
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Fig. 2. Map of integrated emission of the PdBI data (top 
left), and signal-to- noise ratio (top right) for the ^^CO 
(J=l-0) line. Same for the 3mm continuum emission (bot- 
tom panels) . The synthesized beam is shown in the bottom 
left inserts. 



The integrated emission detected with the PdBI is dis- 
played in Fig. [5] (left panel), with the corresponding signal- 
to-noise ratio (right panel). The integrated emission covers 
most of the mosaic area. This is no longer true when this 



emission is displayed in velocity slices (Fig. [31 top panels). 
Several distinct structures are detected in addition to the 
bright CO peak, at velocities [ -3.1, -2.3 ] km s^^. Most are 
weak (the first level in the PdBI channel maps of Fig. [3] 
is 3(t) but they extend over many contiguous synthesized 
beams (10 to 30). 

The PdBI data merged with the short-spacings pro- 
vided by the 30m telescope and the ^■^CO(l-O) emission 
detected by the IRAM-30m telescope are displayed in the 
same velocity-slices, for comparison, in Fig. [31 central and 
bottom panels respectively. Most of the structures seen by 
the PdBI lie at the edge in space and in velocity space of 
extended emission present in the single-dish channel maps. 
This property is most visible for the two structures in the 
north-west of the mosaic over [-4.8, -4.4] km s~^ and [-2, 
-1.2] km s~^, and in the central region &i v — —2.8 km s^-'^. 
It is even better seen by comparing the single-dish maps be- 
fore and after combination with the PdBI data. The single- 
dish maps are changed in two-ways: the structures exhibit 
sharper, more coherent boundaries and these boundaries 
extend further in velocity-space (e.g. channels -4.7 and -2.3 
km s~^). In a given channel of width Awc, the size of the 
detected structures in the CO emission Axc is inversely pro- 
portional to the velocity-shear, Axc — Awc/ (dvLSR./ dxpos)- 
Hence, the detection of small-scale structures at the edge 
of the velocity coverage of larger-scale structures may be 
favored by an increase of the velocity shear at these edges. 

The fact that these structures appear both in PdBI-only 
data and in combined (PdBI+30m) data gives confidence in 
their reality, independently of the deconvolution techniques. 

In summary, the interferometer is sensitive by construc- 
tion to small-scale (i.e. sharp) variations of the space- 
velocity CO distribution. It happens that the sharp struc- 
tures detected by the interferometer lie at the edge in space 
and velocity of regions of shallow CO emission that extend 
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over at least arcminutes, as displayed in the 30m channel 
maps. The PdBI-structures are therefore the sharp edges of 
extended structures. 

3.2. Observed characteristics of the PdBI structures 

We have identified eight structures in the space-velocity 
^^CO(J=1~0) PdBI data cube that are well separated from 
one another in direction and in velocity. They are shown 
in Fig. m each drawn over its proper velocity range. The 
right panels show the PdBI data combined with 30m data 
over the same velocity ranges to further illustrate that the 
PdBI filtering emphasizes the sharpness of the edge of the 
space- velocity structures. Fig. H] also shows that the single- 
dish structures cover a large fraction of the mosaic area. 
For instance, in the case of structure #1, the single-dish 
structure extends over the whole southern half of the mo- 
saic, while for structure #5 it almost covers the northern 
half. 

The observed properties of the 8 PdBI structures are 
given in Table 2. The peak ^^CO(J=1--0) temperature is 
that detected by the PdBI, therefore the excess above the 
extended background, resolved out by the PdBI. The size 
01/2 is the half-power thickness of the elongated structures, 
deconvolved from the beam size. The projected thickness, 
in mpc, is called l± by opposition to the unknown depth 
along the line-of-sight (/os), called /||. The position-angle 
PA is that of the direction defined, within ±10°, by the 
three brightest pixels of each structure. When they are not 
aligned, as in the case of #8, we determine a direction with 
the meaning of a least-square fit. It corresponds to an av- 
erage PA over the detected structure that does not take 
into account the substructure visible in Fig. [S] for instance. 
Because of their different velocity width and CO line tem- 
perature, the CO integrated brightness of the eight struc- 
tures varies by a factor 25. 

Most of the PdBI-structures are elongated and straight 
with different position-angles in the sky. Interestingly, they 
do not shadow each other in space and in velocity space 
(i.e. each fills only a small area of the mosaic in a small 
velocity interval, and the positions and areas of the detected 
structures are different). Their cumulative surface filling 
factor in the mosaic field is large, fs ~ 0.5 (Fig. i.e. 
fs = 0.6 for the structures detected at more than 1-sigma 
and fs = 0.3 for 3-sigma detections. However, the fraction 
of the single-dish power (integrated over the mosaic) seen 
by the PdBI in the CO (1-0) line is low. It depends on the 
velocity interval: it varies between 2% in the ^^CO line-core 
(defined as the velocity range, [-5.0, -3.5]km s^"'^, over which 
the single-dish "CO/^^co is the largest, see F98), and 6% 
in the line- wings. Fig. O displays the emission profile of the 8 
PdBI-structures with the single-dish ^^CO and ^^C0{J=1- 
0) emissions over the same area (defined by the polygons 
of FigHl). 

Last, the PdBI-structures cover the full velocity range 
of the single-dish CO line (see bottom panel of Fig. [5]) in- 
cluding the far line-wings (e.g. structure #2 at -5.5 km 
s~^). Note however that the spectrum integrated over the 
whole mosaic peaks at -3 km s~^, in the wing of the single- 
dish ^^CO line while its minimum, around -4.5 km s~^, 
coincides with the peak of the single-dish ^'^CO line (i.e. 
line core). The broad velocity distribution of the PdBI- 
structures within the single-dish line coverage ensures that 
they are not artefacts of radiative transfer. If they were. 
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overplotted (resp. black and red contours) to display their 
close spatial correlation. 

they would appear preferentially at extreme velocities be- 
cause CO photons escape probability is larger there. There 
may be a small effect since the power fraction in the line- 
wings is slightly larger than in the line-core, but these frac- 
tions are a few percent in each case. The structures found 
are therefore real edges in space and velocity-space of larger 
structures. 

In this respect, it is interesting to place each PdBI- 
structure in its CO (1-0) larger-scale environment at the 
appropriate velocity (Fig. [8]). The PdBI-structures, marked 
as polygons, lie at the edge of structures that extend be- 
yond the field of the mosaic, up to ^ 300 " or 0.2 pc. In 
the case of structures #3, #4 and #5, the orientation of 
the edges of the large-scale patterns is more visible in the 
"CO(l-O) maps (Fig. O, likely because of the 1^00(1-0) 
optical depth. This coincidence strongly suggests that the 
orientation of the PdBI-structures is not only real but also 
rooted in the larger-scale environment. 

3.3. Pairs of parallel structures 

One of the most challenging finding of this study is the fact 
that among the eight elongated PdBI-structures, six form 

3 close pairs (separated by less than 20 " in projection) of 
structures parallel within ±10° (Table 2 and Fig.g]). These 
are the pairs of structures [#3, #8], [#1, #5] and [#6, #7]. 
The average position-angles of each pair PA — 60, 100 and 
168° are all different. Since the structures (at least in the 
two first pairs) are at different velocities, they are not due 
to artefacts of the deconvolution process. 

The probability of a chance association of these three 
pairs in the field of the mosaic is estimated to be at most 

4 X 10~^. It is the cube of the probability of having one close 
pair of parallel structures. The latter is the product of the 
probability, equal to 5.4 x 10~^, that two, out of eight, ran- 
domly oriented straight structures be aligned within ±10° 
of each other (i.e. be together in a solid angle Af2 = 0.1 
sr), by that (ranging between 0.2 and 0.3 depending on the 
orientation of the pair) to be separated in projection by less 
than 20 " in a mosaic of 1' x 2'. The probability of a chance 
association is only slightly underestimated if one considers 
the structure #8 that is not straight, strictly speaking. 



Falgarone et al.: Extreme velocity-shears and CO on milliparsec scale 



60 



C3 
>6 



-50 



Struclus-s j)i 
-5.7 < t < -3.2 km/s 



J I I I I I L 




0.6 




5 



structure J/S 
< y -4-9 km/s 




Struclure #S 
-3.4 < V < -<i.6 km/a 




-eo 



T — I — I — I — ] — I — I — I — r 



Structure 4? 
-3.2 < V < -3 km/n 




5D D 



structure (/E 
-5.6 < » <: -S.4 km/si 




Structure Iti 
-4 .a. < T e. -4,1 km/a 
— ' — ' — ' — ' — i — ' — ' — ' — ' — 




Struelure #6 
-3,4 < v < -<i.6 km/a 




&0 



-] — I — I — I — I — ; — I — I — I — r 



Struclurc 1/6 
-l,f < T < - 1 .3 km/a 
0.6 I ~ 




0.4 



4< 



I 



0.S 



1-S 



0.5 



I 



6t> 



-GO 



50 D 



-50 



Fig. 4. The 8 structures described in Table 1: Left panels: PdBI-only ^^00(1-0) emission integrated over the indicated 
velocity interval appropriate to each structure. Right panels: Same for the combined PdBI+30m mission. The polygons 
show the area over which the CO spectra of FigO are computed. 



Since the probability of a chance association of the ob- 
served pairs is so low, we infer that the pairs are real asso- 
ciations. This physical connexion is supported by the detail 
of the spatial distribution of the ^^CO emission integrated 
over the two velocity ranges of structures #3 and #8 in 
Fig. [5| the hole, in the low-velocity emission is filled in by 
high- velocity emission, while a common average orientation 
exists over ^ 1' for the pair. 

Two position-velocity cuts (Fig. [7|) across the pair 
[#3,#8] further illustrate what is meant by sharp edges and 
real association. The cut across the PdBI-only data cube 
(left panel) shows two CO peaks centered at offset positions 
40 " (resp. 46 ") and velocities -4.9 km s~^ (resp. -1.6 km 
s""'^) for the low- and high- velocity component respectively. 



These resolved peaks are located exactly at the terminal 
pixels of the larger-scale structures visible at the same ve- 
locities in the cut across the PdBI-|-30m data cube (right 
panel). On this cut, the low- velocity component may be fol- 
lowed over all offsets below « 46 ", while the high- velocity 
component is visible at all offsets above w 40 ". This cut 
also illustrates a clear difference between the two velocity 
components: the velocity of structure #3 (peak at -5 km 
s~^ in the PdBI spectrum of Fig. [5]) falls within the veloc- 
ity coverage of the bright extended gas (^^CO and ^'^CO 
line core in the single-dish spectra) while that of structure 
#8 (peak at -1.5 km s^^ in the PdBI spectrum of Fig. [5]) is 
not blended with any other emission in that extreme veloc- 
ity range and appears as a weak emission in the single-dish 
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of 10. Note that i) only a small fraction of the singlc-dish flux is recovered at PdBI and ii) the centroid velocities of the 
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spectrum (i.e. a line- wing). Such blendings in space and ve- 
locity projections with extended components resolved out 
by the PdBI observations (Fig. [71 right panel) explain why 
such pairs of structures are so difficult to recognize in single- 
dish observations or low sensitivity interfcromctric observa- 
tions. 

The PdBI-structures cannot therefore be understood as 
isolated entities. Not only are they the sharp edges of larger 
CO-structures seen in the single-dish maps but also 6 out 
of 8 of these edges are paired. In the following, we will call 



the CO extended structures bounded by sharp edges either 
Sharp-Edged Extended Structures (SEES) or Sharp-Edged 
Extended Double Structures (SEEDS) when they belong to 
a pair, to emphasize this essential property. 

3.4. Velocity shears 

The pairs being real associations, we ascribe a velocity- 
shear to each of them. The projected separation Sl± and ve- 
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locity difference Svlsr between the low- and high- velocity 
components of each pair provide a measure of the velocity- 
shear 6vlsr/^1±- We cannot determine whether this mea- 
sure is a lower or upper limit of the true velocity-shears 
because of the projection effects: both the separation mea- 
sured in the pos and the velocity difference are lower limits. 

The results are given in Table 4. The method used is il- 
lustrated in Fig. [7] (left panel) for the pair [#3,#8]: the 
projected separation between the low- and high-velocity 
components is 6 " or 4.5 mpc while the velocity separation 
is 3.5 km s~^, hence a velocity-shear of 777 km s~^ pc~^, 
the largest ever measured in CO emission in a molecular 
cloud devoided of star formation activity. These values cor- 
respond to an average over several positions along the shear 
direction, including those where the two velocity compo- 
nents partially overlap. The separation is therefore slightly 
underestimated by the averaging. Note that one pair only, 
[#6, #7], has a very small velocity-shear, probably because, 
in that case, the two velocity components involved in the 
shear are mostly in the pos. A rate-of-strain, defined as 
a — ^Svi^sji/6l±, and timescale r = a~^, are also given 
to help comparison with chemical models (Section 7). The 
large observed velocity-shears translate into timescales as 
short as a few lO'^ yr, if the Lagrangian and Eulcrian views 
of the fluid can be exchanged (see .Mordant et al...200^ . 

3.5. The SEE(D)S are layers of CO emission 

The small-scale structures detected by the PdBI have prop- 
erties never seen before because the present observations 
are most sensitive and the field of view is large in compar- 
ison to the resolution: (1) they are not clumps, but elon- 
gated structures, only bounded by the hmited size of the 
mosaic, (2) they all mark a sharp fall-off of the CO emission 
in selected velocity ranges: they are not isolated filaments, 
but the sharp edges (3 to 11 mpc in projection), simulta- 
neously in space and velocity-space, of larger structures, 
the SEE(D)S, extending beyond the mosaic {I > 0.2 pc), 
(3) six of these form three pairs of parallel structures at 
different velocities, with a small projected separation and 
the velocity-shears estimated for two of these pairs, several 
100 km s^^ pc~^, are the largest ever measured in non-star 
forming clouds. 



If the SEE(D)S were CO-emitting volumes (i.e. 3- 
dimensional structures in space) of characteristic dimension 
Z, their edges would be surfaces commensurate with P. In 
projection, these edges would appear as surfaces, also com- 
mensurate with for a random viewing angle. Only if these 
surfaces were plane and viewed edge-on (within ±5 deg for 
a projected size less than one tenth of their real size) would 
these edges appear as thin elongated structures. We rule 
this out on statistical grounds: the mere fact that we de- 
tect 8 sharp CO-edges in the small field-of-view of the PdBI 
observations suggests that it is not a rare configuration and 
that the eight sharp CO-edges are seen from random view- 
ing angles. We thus infer that the SEE(D)S are CO-layers, 
rather than volumes and that their thickness is 10 mpc 
or less, on the order of the width of the PdBI-structures. 

This statistical argument is reinforced by the presence 
of pairs. The SEEDS are structures that have sharp edges 
with only small or null overlaps. If their interface were 2- 
dimensional (i.e. if the SEEDS were volumes), the small 
overlap would occur only for an edge-on viewing, an un- 
likely case. Their interface is therefore 1-dimensional rather 
than 2-dimensional and the SEEDS are layers of CO emis- 
sion. This ensures that under any viewing angle the two 
extended velocity components are detected with only a nar- 
row or null spatial overlap in projection. The SEEDS could 
still be be 3-dimensional pure velocity-structures, where 
large velocity-shears produce sharp edges in channel maps 
of finite spectral resolution (see Section 3.1). However, with 
the same statistical argument as above, concerning density 
structures, we rule out the possibility that the SEEDS be 
5-dimensional velocity structures. These must be CO lay- 
ers. 

In summary, the sharpness of the edges of the SEE(D)S, 
associated with the fact that we detect 8 cases in the mosaic 
and three close-pairs that do not overlap in space, implies 
that the SEE(D)S are thin layers of CO emission rather 
than volumes. 

4. Gas density of the PdBI-structures 

4.1. Estimates from CO line emission 

Because of the elongated shape of most of the structures 
and the fact that they are edges of more extended emis- 
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Table 3. Spatial and kinematic characteristics of the three pairs of parallel PdBI-only structures. 
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length over which the structures are parallel within ±10°. 



sion, we have not tried to decompose the observed emission 
using clump finding algo rithms such as GAUSSCLUMP 
(jStutzki fc Guest"enll990t ). We estimate below the gas den- 
sity in these structures in two independent ways and com- 
pare the results to those inferred from the dust continuum 
emission. 

First, we compute upper limits of the H2 densities 
(Table 2) by adopting the C0-to-H2 conversion factor 
X = 1.56 ±0.05 x 10^0 cm-2 (Kkms"^)"^ (jHunter et all 
119971 ) so that tih^ = 5 x lO"* cm'^ W{CO)/lmpc, for a los 
depth equal to the projected thickness Since we are 
observing edges of layers (see section 3.2), the inferred den- 
sities are overestimated by the unknown factor The 
upper limits of the H2 densities derived from the galactic 
CO to H2 conversion factor (Table 2) vary by a factor 10 
only. 

Alternatively, one may use a LVG analysis to estimate 
the gas properties in these structures. Two assumptions are 
made: (1) the CO emission is not beam-diluted and (2) the 
excitation is assumed to be the same as measured in the 
same field with the IRAM-30m so that we adopt the line 
ratio R(2 - 1/1 - 0) = 0.7± 0.1 (|Falgarone et al.lll998t 
iHilv-Blant et al.ll2008f) . This value may be representative 
of the excitation of translucent molecular gas because the 
same line ratio is found in differe nt observations s ampling 
a similar kind of molecular gas (|Petv et al.l 120081 ). Under 
these conditions, the CO column densities per unit velocity 
are very well determined for all line temperatures. They 
are given in Table 4 for the brightest, weakest, and most 
common CO peak temperature observed. The inferred CO 
column densities differ by only a factor 10 to 16 between 
the brightest and weakest structure. 

Table 4 also gives the range of gas kinetic temperatures 
and associated range of densities, thermal pressures Pth/k 
and CO abundances of possible solutions. The range of tem- 
peratures is bounded towards high values by the thermal 
width of the CO lines {Tk < 250 K for the broadest line, 
< 35K for the narrowest). Solutions colder than lOK are 
unlikely because the gas is poorly shielded from the ambi- 
ent ISRF. The CO optical depth is therefore smaller than 
a few, in agreement with the re sults of Paper II. A similar 
conclusion has been derived bv iHeithausenI (j2006D after he 
failed to detect the ^^CO and Ci8O(J=l-0) line with the 
PdBI in a nearby small-area molecular structure (SAMS) 
field. 

Each set of II2 density and kinetic temperature in the 
LVG solutions, corresponds to a product X(CO)Z|| where 
X{CO) is the CO abundance relative to H2. The range of 
CO abundances inferred from the LVG computations are 



given in Table 2 for ^|| = They may be overestimated 
by the unknown ratio 

Table 4 shows the range of possible II2 densities derived 
from the LVG analysis for gas temperatures between lOK 
and 200K at most. The comparison of these values with the 
upper limits inferred from the CO-to-Il2 conversion factor 
(Table 2) provides narrower II2 density ranges, 71H2 — 800 
to 10^ cm~^, and 300 to 2 x 10^ cm~^, for the weakest 
and most common structures respectively. In spite of all 
the uncertainties, the two methods infer consistent H2 den- 
sities that do not exceed 3 x lO^cm"'^. Moreover, whether 
the gas is cold or warm, its thermal pressure is about the 
same, within a factor of a few, and is in harmony with that 
inferred from carbon line observations in the local ISM that 
has an average oiPth/ fc ~ 3 x 10'^ K cm~^ wit h fiuctuations 
up to ~ 10^ K cm-3 ([Jenkins fc TripDll2007[ ). 

4.2. Estimates from the dust continuum emission 

In addition to ^^CO lines, we have detected continuum 
emission. This emission is close to the noise level, except 
for the large bright spot associated to the ^^CO peak of 
emission. A comparison of the continuum emission with 
the CO contour levels (Fig. ^ suggests that the elongated 
feature of continuum emission in the north-western corner 
of the mosaic is also real. 

On the basis of the coincidence of the peaks of the ^^CO 
and continuum emission in the mosaic, we ascribe the con- 
tinuum emission to thermal dust emission. The average 
continuum brightness over the ^^CO peak (Fig. [5]) is 2±1 
mJy/beam, hence Icont = 1±0.5 x 10"^" ergcm"^ s~^ sr~^ 
Hz"^. The dust opacity, = Icont/ B^iTd), depends on the 
dust temperature. We adopt a dust emissivity = 8.7 x 
10"2S(A/250/im)-2 Nh-, deduced from COBE data for dust 
heated by the ambie nt interstellar radiation field (ISRF) 
()Lagache et al.lll999D . For = 10 K and v = 115GHz, we 
find Td = 3.2 X 10"'^ and Nh = 4.0 x 10^° cm'^ across the 
peak. For = 17 K, the most plausible value in translucent 
gas of the Solar Neighborhood, this value would be lower 
by a factor 2 and for = 8 K it would be larger by 40%, 
providing the range Nh = 2.0 x 10^° to 5.6 x 10^° cm~^ for 
high and low dust temperatures, respectively. This estimate 
of Nh may be compared to that inferred from ^^CO in the 
previous section. If we allow for a column density of atomic 
hydrogen comparable to that of H 2 , as fo und in the Polaris 
Flare bv iHeithausen fc Thaddeui (|l990l ). the total H col- 
umn density inferred from CO for the peak of structure #5 
ranges between iVij = 37V(H2) = 2x10^^ 10 2.4x10^° cm-^ 
for the warm and cold solutions respectively (assuming a 
size of 9mpc) . The two ranges of values would overlap for a 



Falgarone et al.: Extreme velocity-shears and CO on milliparsec scale 11 



Table 4. Results of LVG radiative transfer calculations for representative observed values (see Table 2) 
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depth a few times larger than the observed projected size, 
allowing for warm and moderate density solutions where 
the dust temperature is lower than that of the gas. Given 
the many uncertainties in the different steps (the question- 
able validity of the CO line analysis, the knowledge of dust 
emissivity and temperature) , the consistency between these 
two independent estimates is encouraging and we are con- 
fident that we have detected the dust thermal emission of 
the brightest small-scale structure and that its H2 density 
is not higher than a few lO^cm"^. 

Aside from the peak, the rest of the structures have col- 
umn densities a few times smaller and their dust continuum 
emission is expected to lie closer to the noise level. In addi- 
tion, the surface filling factor of the ^^CO structures being 
large in the central area of the mosaic, the PdBI visibility 
of the continuum emission of individual structures is ex- 
pected to be highly reduced compared to that of the line 
which takes advantage of velocity-space. This may be the 
reason that the continuum emission and the ^^CO emission 
do not coincide elsewhere: the continuum emission is more 
heavily filtered out by the interferometer than the ^^CO 
emission. 

5. What are the SEE(D)S? 

5.1. Manifestations of the small-scale intermittency of 
turbulence 

The two largest velocity- shears given in Table 3 are more 
than two orders of magnitude larger (within the uncertain- 
ties due to projections) than the average value of 1 km 
s~^ pc~^ estimated on the pa rsec scale in molecular clouds 
(|Goldsmith fc Arguillal [l985l ). The velocity field in these 
two SEEDS therefore significantly departs from predictions 
based on scaling laws obtained from CO (1-0) in molecular 
clouds, such as that shown in FiglTOl In spite of a signif- 
icant scatter of the data points, a power law Svi oc l^^^ 
characterizes the increase of the velocity fiuctuations with 
the size-scale at least above ~ 1 pc. Below that scale- 
length, the scatter increases and a slope 1/3 would not be 
inconsistent with the data. According to the former scaling, 
the velocity-shear should increase as l^^^^, therefore by no 
more than 140^^^ = 12 between 1 pc and 7 mpc. If the 
other scaling is adopted, this factor becomes 140^/'^ = 26. 
Now, the observed shears increase by more than two orders 
of magnitude between these two scales. This is conspicu- 
ous on Fig[TO] where the 8 PdBI-structures of Table 2 are 
plotted either individually or as pairs (i.e. as they would be 
characterized if the spatial resolution were poorer and indi- 
vidual structures were not isolated in space, providing for 
instance a linewidth Avi^2 = 3.5 km s^^ and a size l± 7 
mpc for the pair [#3,#8]). 



This result has to be put in the broader perspective 
described in Section 1. The statistical analysis of the ve- 
locity field of this high latitude cloud (Paper III, HF09) 
shows that the pdf of the ^^CO line-centroid velocity in- 
crements increasingly departs from Gaussian as the lags 
over which the increments are measured decrease. The locus 
of the positions that populate the pdf non-Gaussian wings 
forms elongated and thin (~ 0.03 pc) structures that have 
a remarkable coherence, up to more than a parsec. HF09 
propose, on this statistical basis, but also because of their 
thermal and chemical properties given in Section 1, that 
these structures trace the intermittency of turbulent dissi- 
pation in the field. The pair of structures [#3,#8] belongs 
to that locus of positions (see their Fig. 3). The extremely 
large velocity-shears measured in that small field are not 
just exceptional values: they have to be understood as a 
manifestation of the small-scale intermittency of interstel- 
lar turbulence, as studied on statistical grounds in a much 
larger field. 

5.2. The emergence of CO- rich gas 

The PdBI-structures mark sharp edges in the ^^CO emis- 
sion. As discussed in Section 3.3 and illustrated in Fig. [71 
the CO emission of space- velocity structures extending over 
arcminutes (the SEE(D)S) drops below the detection level 
over 4.3 " (the resolution). Therefore, several questions 
arise: what is the nature of the undetected gas that provides 
the continuity of the fiow? Is it undetected because its den- 
sity is too low to excite the J = 1 — transition of ^^CO? 
Or is it dense enough but with too low a CO abundance? 
For simplicity, in the following discussion, "CO-rich" quali- 
fies the gas with X{CO) > 10^^, the CO abundance of the 
weakest detected structure (Table 4), and "CO-poor" the 
gas with a lower CO abundance. 

Acording to LVG calculations, the 3a detection limit of 
our observations allows us to detect CO column densities 
as low as iV(CO) ~ a few 10^^ cm~^ in gas as diluted 
as nn ~ 50cm~"^, at any temperature, and for a velocity 
dispersion of 0.2 km s~^, characteristic of the structures 
found. This detection limit is very low. Therefore, if the 
undetected gas on the other side of the edge is CO-rich 
(with a total hydrogen column density comparable to that 
of the detected part), it has to be at a density lower than 
nn ~ 50cm~'^, not to excite the ^^CO(J=1~0) transition at 
a detectable level. We rule out this possibility because this 
density is that of the CNM and it is unlikely that gas at 
that d ensity be CO-rich (see also the models of iPetv et al.l 
(l200l ). 

The alternative is that the undetected gas is CO-poor 
and that it is not its low density but its low CO abundance 
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Fig. 8. Integrated maps of the same velocity range as that 
defined in Fig.glof the ^^CO (J=l-0) map observed at the 
IRAM-30m. Green crosses delimit the mosaic position and 
the red polygon as defined on Fig. |4] shows the position of 
the elongated structures detected at PdBI. 



that makes it escape detection in ^^CO(J=1-0). Given the 
sharpness of the edges of the SEE(D)S, between 3 and 11 
mpc (Table 2), the process responsible for this transition 
has to be able to generate a significant CO enrichment over 
that small scale. 

In the above, we rule out the possibility that the sharp 
edges (i.e. the PdBI-structures) mark photodissociation 
fronts, because the orientations of such fronts would not 
be randomly distributed, as is observed. Moreover, there 
is no source of UV photons in that high-latitude cloud 
and the radiation field there is the ambient galactic ISRF. 
Photodissociation fronts would not have different orienta- 
tions depending on gas velocities varying by only a few km 
s~^. The sharp edges are not either folds in layers of CO 
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Fig. 9. Same as Fig. |S] except that the single dish map is 
the one of ^^CO (J=l-0). 



emission because those who belong to SEES (single struc- 
tures) lack the second part of the layer, and those who 
belong to SEEDS have the two parts at different velocities. 

We thus infer that the SEE(D)S are the outcome of a 
dynamical process, that involves large velocity-shears, and 
takes place in a gas undetected in CO (1-0) emission, be- 
cause it is CO-poor, not because it is too diluted. This gas 
may be the CNM and the dynamical process has to be able 
to enrich the CNM in CO molecules within a few 10'^ yr 
and over a few milliparsec. 



6. Comparison with other data sets 

Our results broaden the perspective regarding the exis- 
tence of small-scale CO structures in molecular clouds. 
Heithausen (2002, 2004, 2006) has found small-area molec- 
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Fig. 10. Size-linewidth relation for a large sample of 
^^CO(l-O) structures (see Appendix B) to which are 
added: the SAMS data (single-dish data from Heithausen 
f2002, 2006) (soli d triangles), PdBI data within SAMS2 
(jHeithausenl I2004D (open triangles)), a polygon that pro- 
vides the range o f value s for the 12 structures of 
ISakamoto &: Sunadal (|2003f) and the eight structures of 
Table 2 (solid squares). The 3 empty squares without error 
bars show where the three pairs of PdBI-structures would 
be if not resolved spatially (i.e. the velocity increment be- 
tween the two structures would appear as a linewidth for 
the pair). Same with the large triangle for the pair of struc- 
tures in SAMS2. The straight lines show the slopes 1/3 and 
1/2 for comparison. 



ular structures (SAMS) that are truly isolated CO fea- 
tures in the high latitude sky. PdBI observations of the 
SAMS (Heithausen 2004) reveal bright sub-structures that 
are all brighter and broader than our PdBI-structures. 
Unfortunately, the emission has been decomposed into 
clumps, a questionable procedure because short-spacings 
have not been combined to the PdBI data and the CLEAN 
procedure tends to create structures on the beam scale. 
The large H2 densities inferred are therefore likely overes- 
timated. An interesting feature can be seen in the channel 
maps, though. Two elongated thin patterns cross the field, 
reminiscent for their thickness and length of what is found 
in the present study. A velocity shear of 180 km pc~^ 
is estimated between these two elongated structures for a 
velocity separation of 0.9 km and a pos spatial separa- 
tion of 10 " on average (or 5 mpc at the assumed distance 
of 100 pc). This velocity-shear is thus commensurable with 
the two largest values found in the Polaris field. 

Ingalls et al. (2007) have detected milliparsec clumps in 
a high latitude cloud. They are located in the line-wings of 
the CO single-dish spectrum and they model them as tiny 
(1-5 mpc) clumps of density of a few lO'^ cm~^. A more 
detailed comparison with the present results is not possible 
because they do not analyze individual structures. 

Sakamoto & Sunada (2003) have discovered a number 
of CO small-scale structures in the low-obscuration regions 
of long strip maps beyond the edge of the Taurus molecular 
cloud. Their main characteristics are their large line-width 
and their sudden appearance, and disappearance, within 
0.03 to 0.1 pc. The authors interpret these features as the 
signature of structure formation induced by the thermal in- 



stability of the warm neutral medium (WNM) in the turbu- 
lent cloud envelope. These CO small-scale structures thus 
carry the kinematic signatures of the embedding WNM, 
hence their large velocity dispersion, both interclump and 
intraclump. The inferred line ratio, i?(2 — 1)/(1 — 0) = 0.4, 
is low, consistent with a low excitation temperature and H2 
densities lower than ~ 10^ cm"'^. The authors propose that 
their small-scale CO structures pinpoint molecule-forming 
regions, driven by the thermal instability in the turbulent 
diffuse ISM. 

Our data therefore share many properties with these 
different samples. FigdO] allows a comparison of the pro- 
jected size and linewidth of the above milliparsec-scale 
structures with those of CO (1-0) structures identified in 
data cubes from non-star-forming regions of all sizes, up to 
several 100 pc (see the relevant references in Appendix B). 
Although some of them (a few individual PdBI-structures 
of our sample) further extend the general scaling law down 
to 2 mpc, most of them significantly depart from this law 
by a large factor. As already mentioned in Section 5.1, the 
departure is the largest for the pairs of PdBI-structures, as 
they would appear if they were not resolved spatially i.e. 
as anomalously broad structures with respect to their size. 
The increased scatter of velocity-widths of the structures 
below 0.1 pc down to 1 mpc may be seen as another man- 
ifestation of the intermittency of turbulence in translucent 
molecular gas. 

7. Comparison with experiments, numerical 
simulations and chemical models 

The present data set discloses small-scale structures of in- 
tense velocity-shears that carry the statistical properties 
of intermittency and, in conjunction with that of IIF09, 
reveals a connexion between parsec-scale and milliparsec 
scale velocity- shears. The dynamic range of coupled scales 
in the Polaris Flare is therefore on the order of ~ 10'^. 
Moreover, velocity differences, up to 3.5 km s^^, close to 
the rms velocity dispersion of the CNM turbulence mea- 
sured on 10-pc scales (or ni ore) (|Miville-Deschenes et al.l 
[2OO3I : iHaud fc Kalbin^l200l . are found in the PdBI field 
over 10 mpc, without any detected density enhancement 
nor shock signature. We argue that the SEE(D)S are the 
CO-rich parts of straining sheets in a gas undetected in 
1^00(1-0), likely the CNM, and that the fast CO enrich- 
ment is driven by enhanced turbulent dissipation in the 
intense velocity-shears. We show below that these findings 
may be understood in the light of recent numerical simu- 
lations of incompressible an d compressible turbu lence, and 
the TDR chemical model of iGodard et~all (|2009t l. 

The fact that the most dissipative structures appear to 
be layers of intense strain-rate is consistent with recent re- 
sults of numerical simulations of i ncompressible turbu lence 
at high Reynolds num ber (iMoisv fc Jimene 3120041) an d 
laboratory experiments (jCanapathisubramani et aLfr2008l ) . 
These regions are not randomly distributed and form 
inertial-range clusters (jMo isv & Jimenez 2004) or develop 
at the bound aries regions of high level of vo rticity (i.e. 
vortex tubes) ([Ganapathisubramani et al.ll2008( l. Coupling 
between small-scale statistics of the velocity field and the 
properties of the large-scale flows is al so clearly pr obed in 
the high-i?e numerical simulations of Mininni et al.l (2006): 
correlations are observed between large-scale shear and 
small-scale intermittency. 
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In compressible turbulence, the fact that the most dis- 
sipative structures are shear-layers is not expected. Yet, 
in their hydrod ynamical sim u lation s of mildly compress- 
ible turbulence, iPorter et alJ (|2002f l show that the com- 
pressible component of the velocity field is weaker than 
its solenoidal counterpart by a factor ^ 3, independent 
of the nature of the drivin g proc ess (compressible or 
solenoidal) and IVestuto et all (|2003[ ) find that the energy 
fraction in the solenoidal modes is dominant and increases 
with the magnetic field intensity in compressible magneto- 
hydrodynamical (MHD) turbulence. These numerical ex- 
periments are still far from approaching the ISM conditions 
but they suggest that turbulent dissipation may occur pri- 
marily in solenoidal modes, i.e. without direct gas compres- 
sion, and that the properties of the small scales are coupled 
to the large-scales. 

In the TDR models of iGodard et all (|2009t) . the chem- 
ical enrichment of the CNM is driven by high gas tem- 
peratures and enhanced ion-neutral drift, without den- 
sity enhancement. The temperature increase is due to vis- 
cous dissipation in the layers of largest velocity-shears at 
the boundaries of coherent vortice^. The large ion-neutral 
drift occurs in the layers of largest rotational velocity in 
which ions and magnetic fields decouple from neutrals. 
These two dissipative processes trigger endothermic chem- 
ical reactions, blocked at the low temperature of the CNM. 
Enrichments consistent with observations are obtained for 
turbulent rates-of-strain a = 10~^^ s~^ induced by large 
scale turbulence and for moderately dense gas (nn < 200 
cm~^) characteristic of the CNM. There is no direct de- 
termination of the rates-of-strain generated by turbulence 
in the CNM. We note however that the largest observed 
velocity-shear (Table 3) corresponds, if the projected quan- 
tities provide reasonable estimates, to a comparable rate- 
of-strain. The range of observed CO column densities from 
A^(CO) = 10^^ to 1.6 X 10^^ cm~^ can be reproduced by in- 
tense velocity- shears occurring in gas of density 100 to 200 
cm~^. In this framework, the energy source tapped to en- 
rich the medium in molecules is the supersonic turbulence 
of the CNM. 

The association between the large observed velocity- 
shears and local enhanced dissipation rate is therefore sup- 
ported not only by the earlier works presented in the 
Introduction but also by a quantitative agreement between 
the TDR chemical models and the present observational 
results. We cannot rule out however a contribution of low- 
velocity C-shocks to the turbulent dissipation. If they prop- 
agate in the CNM, they are not visible in the CO lines. 
Such shocks are not yet reliably modelled (Hily-Blant et al. 
in preparation). 

8. Conclusions and perspectives 

IRAM-PdBI observations of a mosaic of 13 fields in the tur- 
bulent environment of a low-mass dense core have disclosed 
small and weak ^^CO(l-O) structures in translucent molec- 
ular gas. They are straight and elongated structures but 
they are not filaments because, once merged with short- 
spacings data, the PdBI- structures appear as the sharp 
edges of larger-scale structures. Their thickness is as small 
as « 3 mpc (600 AU), and their length, up to 70 mpc, 

^ the "sinews of turbulence" put forward by iMoffatt et al.l 

ij 19941 ) that link large-scale strain and small-scale vorticity 



is only limited by the size of the mosaic. Their CO col- 
umn density is a well determined quantity for the excita- 
tion conditions found at larger scale and is in the range 
iV(CO) = lO^-* to 10^5 cm-2. Their H2 density, estimated 
in several ways, including the continuum emission of the 
brightest structure, does not exceed a few 10'^ cm~'^. Their 
well-distributed orientations can be followed in the larger- 
scale environnement of the field. Six of them form three 
pairs of quasi-parallel structures, physically related. The 
velocity-shears estimated for the three pairs include the 
largest ever measured in non-star-forming clouds (up to 780 
km s^^ pc^^). 

The PdBI-structures are therefore not isolated and are 
the edges of so-called SEE(D)S for sharp-edged extended 
(double) structures. We show that the SEE(D)S are thin 
layers of CO-rich gas and that their sharp edges pinpoint a 
small-scale dynamical process, at the origin of the CO con- 
trast detected by the PdBI. We propose that the SEE(D)S 
are the outcomes of the chemical enrichment driven by in- 
tense dissipation occurring in large velocity-shears and that 
they are CO-rich layers swept along by the straining field 
of CNM turbulence. 

The present work is the first detection of mpc-scale in- 
tense velocity- shears belonging to a parsec-scale shear. The 
large departure from average of the kinematic properties 
of these structures, confirms that they are a manifestation 
of the small-scale intermittency of turbulence in this high 
latitude field, a property already established on statisti- 
cal grounds (IIF09). The values of the velocity-shears (or 
rate-of-strain) provide a quantitative constraint on the dis- 
sipation rate that can be compared to chemical models. 
The link between the turbulent dissipation in the diffuse 
gas and the dense core observed in the vicinity of the PdBI 
mosaic (Fig. [T]) still remains to be established. 

Last, we would like to stress that sub-structure still ex- 
ists in these mpc-scale structures of the diffuse ISM and 
that the next generation of interferometers (e.g. ALMA) 
should be able to observe gas at the dissipation scale of 
turbulence (that is still unknown) or at least observe the ef- 
fects on the ISM (temperature, excitation, molecular abun- 
dances) of the huge release of energy expected to occur 
there. 
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Fig. A.l. Map of the noise level in Kkms over the 13- 
field mosaic. 
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Appendix A: Noise level in the mosaic 

Mosaic noise is inhomogeneous due to primary beam cor- 
rection. This is shown in Fig lA.ll The 13-field mosaic pro- 
duces a large area with uniform noise level. Only at the 
edge of the mosaic does it increase sharply due to the pri- 
mary beam correction (the contour shown are at a 2, 3 and 
4 sigma level, 1 sigma being measured at the map center 
on a channel devoided of signal). 



Appendix B: The size-linewidth scaling law 

Molecular cloud parameters have long been determined 
as those of three-dimensional structures isolated in the 
four-dimensional space of the molecular line data sets 
Ti(a;,y,fz), the line brightness temperature being a func- 
tion of position in the pos (two coordinates x,y), and 
one spectral dimension, the projected velocity on the los 
direction Vz- In this 4D space, 3D structure s are iso- 
lated following differe n t methods ptutzki fc Gue sten 199i 



'Williams et al.' '1994"; 'Falgar one fc Peraultl 1 1987c ,Lore" 



1989; Falgarone et al. 1992). The size and linewidth of 
the large number of clouds displayed in FiglTO] have been 
obtained by using published values, corrected in several 
cases to allow the size and li newidth obey th e same def- 
initions in all the samples (see lFalgaroii3ll998f ) . The struc- 
tures are identified in ^^CO(l-O) molecular hne surveys 
of the central parts of the Galaxy (stars , iDame et alj 
(|1986D . open triangles ISolomon et al.l ()1987[ )') and of the 
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third quadrant (open hexagons. iMav et al.l (|1997f )). in the 
Rosette (crosses) and Maddalena (open squares) molecular 
clouds ( jWiUianis et al., ,1994), in non-star- forming clouds 
(solid triangles. iFalgarone fc Peraiiitl (119871), so l id squ ares, 
iFalgarone et all (|1992D . tripo ds. iLenime et all (|1995D ). in 
p Ophiuchus (soHd hexagons. iLqren (I1989D) and in a high 
latitude cloud (starred triangles. iHeithausen et al.] (|l998f ii). 



